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(57) Abstract* Method and apparatus f r detecting and locating leaks in a fluid carrying pipe using correlation based techniques 
is described. Two input signals are detected from the fluid carrying pipe, the phase f the input signals in the frequency domain is 
analysed to provide at least one filter; the input signals are filtered and cross-correlation is performed n the filtered signals. 
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CORRELATION ANALYSIS IN THE PHASE DOMAIN. 

This invention relates to the detection of a common signal within two (or more) 
input signals using correlation based techniques, and is more specifically, 
although not exclusively, concerned with the detection and location of leaks in 
water mains using correlation based techniques. 

There are many applications in which signals are analysed using the technique 
of crosscorrelation of two data streams. This technique is specifically useful 
when analysing or comparing two or more composite time sequential signals 
within which it is believed that there is a common signal. 

One area in which crosscorrelation is used is in the analysis of water leak 
sounds conducted along water mains, specifically locating leaks on water 
mains using audio sensors. The technique involves the detection and location of 
a common leak sound in the signals from two listening devices. It is anticipated 
that the sound of a leak heard remotely at one location will contain significant 
similarities to the sound of the same leak heard remotely at a second location. 
The two audio signals are therefore likely to exhibit a peak in their 
crosscorrelation function if there is a common source present, the position of 
the peak giving the time delay when the signals are similar. The peaks in the 
crosscorrelation function can therefore be used to deduce the location of 
common audio signals via the velocity of the sound in the pipe and the physical 
surroundings of the noise source. If the nature and persistency of the sound is 
indicative of leaks then there is a high probability that a leak has been found 
and located. 
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The analysis techniques described herein both for the prior art methods and the 
present invention are also applicable in many other technical areas. These 
include the following: the analysis of odd mechanical noise in equipment; the 
assessment of quality of objects by the noise that they make while travelling 
along a production line; the assessment of the performance of wheeled vehicles 
on road or railway tracks; the analysis of road or rail surfaces by the analysis of 
the sound of the passage of vehicles either on a mobile platform or at a fixed 
location; and the measurement of closely similar recordings (the analysis of 
their true relationship from a high quality digital copy to a different 
performance of the same song that sounds very similar). 

Correlation processing techniques are frequently performed using the following 
process: calculate the Fourier transform of the data sets representing the two 
waveforms to be correlated, thereby transforrning the data into the frequency 
domain; performing the required multiplications; and calculating the inverse 
Fourier transform to display the correlation function in the time domain. Often 
the Fast Fourier Transform is used since this allows a considerable reduction in 
the processing required to generate the correlation function. These techniques 
are well known and generally used. 

The general correlation processing techniques work well, and several devices 
are available to perform this task. However, when the signal to noise ratio is 
poor and/or the bandwidth of the signals is restricted due to the propagation of 
the signal then the crosscorrelation function inevitably becomes degraded, with 
the peaks becoming less distinct or even hidden by the noise. Thus filters are 
used to obtain the optimum signal to noise ratio and correlation accuracy. 
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Noise which may be present in a signal can be considered as two distinct 
forms: random noise and clutter. Random noise is unpredictable on a short 
timescale, but has a zero mean. It may be removed by averaging over a 
sufficiently long period. Clutter is a function of the detected signal such as, for 
example, echoes and resonance, which cannot be removed by averaging. 

The signal to noise ratio may be improved by filtering the input signals to 
reduce the bandwidth. However, as the bandwidth is reduced the width of the 
correlation peak tends to increase. This makes it difficult to determine the time 
delay when the signals are similar, and thereby makes it more difficult to 
determine the location of the common signal. It is therefore important to filter 
the signal to exclude only those bands which do not contain the useful signal. 
However, the prior art methods of selecting the bands are imprecise. 

For example, in the case of water leak detection the audio noise signal picked 
up from the pipe will be a combination of the leak noise and other external or 
waterborne sounds. In order to improve the signal to noise ratio of the input 
signals the operator of a water leak detector may analyse the spectrum of the 
input signals to determine where signal power is significant. Filters are then set 
to remove the frequencies in which the signal power is not significant. This 
method usually provide some improvements, but does not guarantee a band 
with good correlation. The operator will generally need to use a combination of 
experience and trial and error to improve the shape of the correlation function. 
This process can be time consuming, often taking tens of minutes at each leak 
site. 
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The present invention seeks to mitigate these disadvantages, and to provide an 
improved method of detection and location of a signal using correlation based 
techniques. 

According to a first aspect of the invention, there is provided a method for 
detecting and locating a common signal within two input signals using 
correlation based techniques, comprising providing at least one filter by 
analysing the phase of the input signals in the frequency domain; filtering the 
input signals in the frequency domain using said at least one filter; and 
performing crosscorrelation of the input signals. 

According to a second aspect of the invention, there is provided a method for 
detecting and locating leaks in a fluid carrying pipe such as, for example, a 
water main using correlation based techniques, comprising detecting two input 
signals from the fluid carrying pipe; analysing the phase of the input signals in 
the frequency domain to provide at least one filter; filtering the input signals in 
the frequency domain using the at least one filter; and performing 
crosscorrelation of the input signals. 

According to a third aspect of the invention, there is provided apparatus 
for detecting and locating a common signal within two input signals using 
correlation based techniques, comprising a computer including: means for 
providing at least one filter by analysing the phase of the input signals in the 
frequency domain; means for filtering the input signals in the frequency 
domain using said at least one filter; and means for performing crosscorrelation 
of the filtered signals. 
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According to a fourth aspect of the invention, there is provided apparatus for 
detecting and locating leaks in a fluid carrying pipe using correlation based 
techniques, comprising: detectors for detecting two input signals from the fluid 
carrying pipe; a computer including means for analysing the phase of the input 
signals in the frequency domain to provide at least one filter; means for 
filtering the input signals in the frequency domain using the at least one filter; 
and means for performing crosscorrelation of the filtered signals. 

The signals may be audio signals. 

The at least one filter may include a first filter for suppressing frequencies 
which do not exhibit a sufficient degree of coherence. The first filter may be 
constructed using a method comprising: selecting at least one section from each 
of the two input signals; calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference between the two 
input signals for each of a plurality of frequencies; and calculating the 
magnitude of the vector sum for each frequency. 

The at least one filter may include a second filter for identifying regions in a 
frequency spectrum of a cross correlation function likely to exhibit a correlated 
phase between adjacent frequencies in its fourier transform. The second filter 
may be constructed using a method comprising: selecting at least one section 
from each of the two input signals; calculating the fourier transform for each 
section; calculating the average vector sum of the phase difference between the 
two input signals for each of a plurality of frequencies; and calculating the 
magnitude of the vector sum for each frequency. 
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The time delay between the common signal in the input signals may be 
calculated by tracking the phase difference between the input signals as a 
functional frequency using the second filter. The variation in the time delay 
between the common signal in the input signals as a functional frequency may 
be calculated using a second filter. 

There may be a third filter for removing frequencies which do not have 
sufficient amplitude. The third filter may be constructed using a method 
comprising applying a digital threshold to the product of the spectra of the two 
input signals. 

The at least one filter may include a fourth filter for compensating the input 
signals for dispersion effects. 

A specific embodiment according to the invention is provided for the detection 
of leaks in water mains. The apparatus may comprise a novel correlator 
provided on a standard field engineer's Personal Computer. The novel 
correlator may include a suite of signal processing algorithms and a suitable 
Graphical User Interface. Like current correlators the novel correlator would 
receive its inputs from two transducers. Input signals are processed by the 
computer using signal processing techniques including those described herein 
to detect and locate leaks. 

As explained hereinbefore, the characteristics of a correlated signal are that it 
has similar features in its waveform (although one may be delayed with respect 
to each other and the whole buried in noise and interference). The 
crosscorrelation of two input signals will only produce a significant peak if the 
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frequency bands which make up those input signals have a stable phase 
relationship with respect to each other, i.e. the input signals have a narrow band 
coherence. Filtering of the two input signals to remove or suppress the 
frequencies which do not exhibit a sufficient degree of coherence will therefore 
enhance the peak in the crosscorrelation function. Frequencies which do not 
exhibit sufficient coherence need not be included in the final result because 
they will only be contributing noise. The filter profile for the signal may 
therefore be determined automatically. 

The novel correlator includes a phase confidence filter, which is constructed to 
suppress the frequencies in the input signals which do not exhibit a sufficient 
degree of coherence. The phase confidence filter effectively indicates the 
estimation of how well a particular frequency will correlate, and is the result of 
analysing the phase stability between the input signals. One possible method 
for constructing the phase confidence filter would be first to calculate a phase 
confidence function using the following steps: to calculate the complex Fourier 
transform for each of a number of sections of the two input signals; to calculate 
the average vector sum of the phase difference between the two input signals 
for each frequency; and then to calculate the magnitude of the vector sum and 
normalise the result. Thus a function is provided which has a value of 1.0 for a 
given frequency which shows perfect coherence, and tends to l/>/n if the 
phase differences for a given frequency are random, where n is the number of 
samples taken. The phase confidence function may then be used as an optimal 
frequency weighting function to construct a filter which will include only those 
frequencies which contribute usefully to the crosscorrelation of the input 
signals. 
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The phase confidence filter is used in combination with an automatic band 
filter. The automatic band filter is constructed based purely on the amplitude of 
the input signals. The filter is designed to select frequencies of large amplitude, 
which are likely to be prominent in the crosscorrelation, and by reducing the 
bandwidth to hence improve the signal to noise ratio. This filter may be 
constructed by simply applying a digital threshold to the product of the spectra 
of the input signals. The threshold may be configurable by the user. The phase 
confidence filter and the automatic band filter are applied separately to the 
data. 

The use of the automatic band filter addresses one practical disadvantage of the 
phase confidence filter. This is that some of the noise in the input signals can 
have highly correlated noise as a result of imperfect digitisation, this effect, 
seen usually at high frequencies, can produce a substantial phase confidence 
measurement. The automatic band filter suppresses frequencies where the 
signal to noise ratio is poor and this artifact is most likely to be seen. 

The effect of using the phase confidence filter in combination with the 
automatic band filter is shown in Figs. 1 to 4. Figs. 1 and 2 show a first input 
signal 1 and a second input signal 2 respectively, both in the frequency domain. 
The crosscorrelation function 7 for the unfiltered first and second signals is 
shown in Fig. 3. The phase confidence filter 3 and the automatic band filter 4 
are constructed for the first and second signals and applied to the signals in the 
frequency domain to produce a filtered first signal 5 and a filtered second 
signal 6. The crosscorrelation function 8 for the filtered first and second signals 
is shown in Fig. 4. The improvement between the crosscorrelation functions 7 
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and 8 can clearly be seen. There is improved clarity, with reduced background 
level and less ambiguous peak location. 

If the common features in the input signals contain some degree of time 
dependency, then the phase relationship will not be stable. The use of the phase 
confidence filter calculated using the method above may blur the results. An 
alternative method for calculating the phase confidence filter would be to 
consider the phases for each individual section of the input signals. The process 
could then adapt if the phase appeared to be stable to lengthen the average by 
including further sections to improve the estimate. 

This method could also be adapted to track the signals if they were time 
dependent. Consider, for example, a signal which is very noisy but has a 
particular value which may be either constant or changing over time. If this 
value is constant then it is possible to average the signal to produce the value 
accurately. However, if the value is changing then a compromise must be made 
between the number of samples over which the signal is averaged, and the rate 
at which the value is varying. 

If to a first approximation this variation is linear it is relatively easy to produce 
a root mean square best straight line fit that continually updates as the data 
comes in. If the variation is non-linear then some means must be found of 
releasing the linear constraint but still producing a good curve fit. If the 
calculation process also produces two alternative straight line fits, one based on 
the first half of the data, and the other based on the second half of the data (they 
could be updated every other data point to ensure an even number of points) if 
the two lines show coherent deviation above some threshold a track deviation 
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can be initiated. Thereby, an adaptive signal tracking process may be 
developed. 

The novel correlator also includes a filter based on a consideration of the phase 
relationships between adjacent frequencies. Relationships exist between 
adjacent frequencies in the different phase analyses. That is, as well as having 
a reasonably stable value, the phase differences may exhibit a steady 
progression of phase over the range of frequencies which are coherent. The 
source of this result can be visualised by considering a single correlation peak, 
a peak in the time domain is created by frequencies adding up in phase in the 
frequency domain. Consider firstly a centre correlation peak (often the result of 
an instrumental artefact or interfering signal) which is located in the centre of 
the correlation range and which corresponds to correlation where there is zero 
relative delay between two input signals. In this case the frequencies have the 
same or very similar phase, and therefore the phase does not change with 
frequency. A centre correlation peak may be removed by constructing a filter to 
remove frequencies having the same phase. Note, however, that exact centre 
correlation in the leak, although unlikely, is not impossible. 

If apeak in the cross correlation function represents a delay between the input 
signals then frequencies in the cross correlation function of those signals are 
shifted in phase by an amount that is proportional to the time delay and the 
frequency. Frequencies evenly spaced in a Fourier Transform will exhibit a 
linear change of phase between adjacent correlated phases. 

Considering the phases of sequential frequencies can have the advantage when 
the propagation bandwidth is very narrow resulting in a broad correlation peak 
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because the phase progression can be "unwrapped" and a straight line fitted to 
it. 

The existence of a second major peak in the correlation result will add a second 
linear progression to the phase partem. The two phase progressions will add as 
vectors with the result that as the relative phases of the two linear progressions 
move in and out of phase the main peak will be advanced and retarded in phase 
giving a ripple in the phase results. 

Alternatively, the case where the velocity of the sound is frequency dependent 
may be utilised. In this case the phase progression for adjacent frequencies in a 
single correlation peak will not be linear. This non-linearity may be tracked. 
This will allow either the compensation of results for dispersion where the 
correlation is spread by its effects or the identification and measurement of the 
dispersion effects using known sound sources. 
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CLAIMS. 

1 . A method for detecting and locating a common signal within two input 
signals using correlation based techniques, comprising providing at least 
one filter by analysing the phase of the input signals in the frequency 
domain; filtering the input signals in the frequency domain using said at 
least one filter; and performing crosscorrelation of the filtered signals. 

2. A method for detecting and locating leaks in a fluid carrying pipe using 
correlation based techniques, comprising: detecting two input signals from 
the fluid carrying pipe; analysing the phase of the input signals in the 
frequency domain to provide at least one filter; filtering the input signals in 
the frequency domain using the at least one filter; and performing 
crosscorrelation of the filtered signals. 

3. A method according to claim 1 or 2, wherein the signals are audio signals. 

4. A method according to any preceding claim, wherein the at least one filter 
includes a first filter for suppressing frequencies which do not exhibit a 
sufficient degree of coherence. 

5. A method according to claim 4, wherein the first filter is constructed using a 
method comprising: selecting at least one section from each of the two input 
signals; calculating the Fourier Transform for each section; calculating the 
average vector sum of the phase difference between the two input signals 
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for each of a plurality of frequencies; and calculating the magnitude of the 
vector sum for each frequency. 

6. A method according to any preceding claim, wherein the at least one filter 
includes a second filter for identifying regions in the frequency spectrum of 
a crosscorrelation function likely to exhibit a correlated phase between 
adjacent frequencies in its Fourier Transform. 

7. A method according to claim 6, wherein the second filter is constructed 
using a method comprising: selecting at least one section from each of the 
two input signals; calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference between the two 
input signals for each of a plurality of frequencies; and calculating the 
magnitude of the vector sum for each frequency. 

8. A method according to claim 6 or 7, including calculating the time delay 
between the common signal in the inout signals by tracking the phase 
difference between the input signals as a function of frequency using the 
second filter. 

9. A method according to any of claims 6 to 8, including calculating variations 
in the time delay between the common signal in the input signals as a 
function of frequency using the second filter. 

10. A method according to any preceding claim, including using a third filter to 
remove frequencies which do not have sufficient amplitude. 
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1 1 . A method according to claim 10, wherein the third filter is constructed 
using a method comprising: applying a digital threshold to the product of 
the spectra of the two input signals. 

12. A method according to any preceding claim, wherein the at least one filter 
includes a fourth filter for compensating the input signals for dispersion 
effects. 

13. Apparatus for detecting and locating a common signal within two input 
signals using correlation based techniques, comprising a computer 
including: means for providing at least one filter by analysing the phase of 
the input signals in the frequency domain; means for filtering the input 
signals in the frequency domain using said at least one filter; and means for 
performing crosscorrelation of the filtered signals. 

14. Apparatus for detecting and locating leaks in a fluid carrying pipe using 
correlation based techniques, comprising: detectors for detecting two input 
signals from the fluid carrying pipe; a computer including means for 
analysing the phase of the input signals in the frequency domain to provide 
at least one filter; means for filtering the input signals in the frequency 
domain using the at least one filter; and means for performing 
crosscorrelation of the filtered signals. 

15. Apparatus according to claim 13 or 14, wherein the signals are audio 
signals. 
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16. Apparatus according to any of claims 13 to 15, wherein the at least one 
filter includes a first filter for suppressing frequencies which do not exhibit 
a sufficient degree of coherence. 

17. Apparatus according to claim 16, wherein the first filter is constructed using 
a method comprising: selecting at least one section from each of the two 
input signals; calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference between the two 
input signals for each of a plurality of frequencies; and calculating the_ 
magnitude of the vector sum for each frequency. 

18. An apparatus according to any of claims 13 to 17, wherein the at least one 
filter includes a second filter for identifying regions in the frequency 
spectrum of a crosscorrelation function likely to exhibit a correlated phase 
between adjacent frequencies in its Fourier Transform. 

19. An apparatus according to claim 1 8, wherein the second filter is constructed 
using a method comprising: selecting at least one section from each of the 
two input signals; calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference between the two 
input signals for each of a plurality of frequencies; and calculating the 
magnitude of the vector sum for each frequency. 

20. An apparatus according to claim 18 or 19, including calculating the time 
delay between the common signal in the inout signals by tracking the phase 
difference between the input signals as a function of frequency using the 
second filter. 
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21. An apparatus according to any of claims 18 to 20, including calculating 
variations in the time delay between the common signal in the input signals 
as a function of frequency using the second filter. 

22. Apparatus according to any of claims 13 to 21, including a third filter to 
remove frequencies which do not have sufficient amplitude. 

23. Apparatus according to claim 22, wherein the third filter is constructed 
using a method comprising: applying a digital threshold to the product of 
the spectra of the two input signals. 

24. Apparatus according to any of claims 13 to 23, wherein the at least one 
filter includes a fourth filter for compensating the input signals for 
dispersion effects. 
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